Comment on: Failure of the work-Hamiltonian connection for free energy calculations 



O 

o 

(N 
C 

3 



A. Imparato^B and L. Peliti 2 '0 

1 Dipartimento di Fisica and CNISM, Politecnico di Torino, 
Corso Duca degli Abruzzi 24, 10121 Torino (Italy) 
INFN, Sezione di Torino, Torino (Italy) 
2 Dipartimento di Scienze Fisiche and CNISM, Universita "Federico II" 
Complesso Monte S. Angelo, 80126 Napoli (Italy) 
INFN, Sezione di Napoli, Napoli (Italy) 
(Dated: June 8, 2007) 

We argue that the apparent failure of the work-Hamiltonian connection for free energy calculations 
reported by Vilar and Rubf (cond-mat arXiv:0704.0761v2) stems from their incorrect expression for 
the work. 
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In a recent preprint, Vilar and Rubi [l| have argued 
that the pivotal connection between the microscopic work 
W performed by a time-dependent force and the energy 
change associated with its change in state cannot be used 
to estimate free-energy changes. Since in several occa- 
sions some groups have exploited the measure of work 
applied on microscopic systems to evaluate its free en- 
ergy landscape, the failure of this connection would have 
far-reaching consequences. 

In the present note we reevaluate Vilar and Rubi's ar- 
guments to test their validity. We find that they rest 
on an incorrect expression of the work to be evaluated, 
and that therefore, if one measures the correct quantity, 
which is experimentally accessible, it is possible to exploit 
the work-free energy relation in order to gather thermo- 
dynamical information on microscopic systems. 

Let us consider a system whose microscopic state is de- 
scribed by the coordinate x, and governed by the hamil- 
tonian H(x,ii) which depend on an external parameter 
fi which can be manipulated. The equilibrium state for a 
given value of \x is described by the canonical distribution 



p7( x ) = 



-H(x,y.)/T 



where the partition function is given by 
Z„ = fdx e~ H ^/ T . 



(1) 



(2) 



The free energy G^ associated with the given value of \i 
is given by 



G,, 



-TlogZ,,. 



(3) 



Let us now evaluate the change in G^ as \x varies from, 
say, /j, = to a final value /i. One has 



(-) 



(4) 



where (A) denotes the average of the function A(x) with 
respect to the canonical distribution (JTJ. Thus 



AG — — Go — 



, /d_H\ 



(5) 



We now define the reversible work W Iev as the quantity 
which appears on the rhs of this expression: 



Wrev = / dWrev = dfl 



,/8H\ 



(6) 



This definition is an immediate consequence of equa- 
tion (128.1), p. 535 of Tolman's book [2|. Its contents 
is explained by Tolman [2j, p. 527] in the following way: 

The concept of work performed by a ther- 
modynamic system on its surroundings de- 
pends on the possibility of changing the val- 
ues of external parameters for the system, 
which have the nature of generalized coor- 
dinates. In the thermodynamic treatment, 
when a small variation is made in the value 
of such an external coordinate, the work done 
can be set equal to the product of that vari- 
ation by the corresponding generalized force 
exerted by the system on its surroundings, 
and can be regarded as equal to the (poten- 
tial) energy thereby transferred to external 
bodies. [. . . ] In the corresponding statistical 
treatment we can assign to all the systems in 
the representative ensemble the same value of 
the external coordinates as those of the ther- 
modynamic system of interest, and can repre- 
sent the corresponding generalized forces by 
their average values in that ensemble. The 
work done and the energy thereby transferred 
to external bodies can then be represented by 
its average value for the ensemble, as given 
by the product of the definite value of the 
displacement with the average value of the 
corresponding generalized force. 
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With this definition, the equality between AG and W n 
is trivially satisfied: 



AG = W„ 



(7) 



It is important to remark that, by definition, the re- 
versible work is a non-fluctuating quantity (the equilib- 
rium average of an observable). Thus, if we evaluate, 



e.g. 



-Wrev/T 



), we have 



e W rev /T) =e -W x ^/T = exp . 



On the other hand, because of 0, one has 



-W Iev /T 



-AG/T 



and, because of ©, 



e -AG/T = e -( GfI -Go)/T = 



(8) 
(9) 

(10) 



Equation © represents the special case of Jarzynski's 
equality [3[ which holds in the limit of reversible ma- 
nipulation. Let us define the fluctuating infinitesimal 
work, dW(x, /i) as the quantity whose average is given 
by dW rev (n): 



dW(x, fi) = dfi 



dH(x, n) 
9/i 



Then the Jarzynsky equality states that 



In this expression, W is the fluctuating total work 



(11) 



(12) 



W = dW = 



t'=0 



(13) 

and the average is taken over all realizations of the pro- 
cess x(t) with a given manipulation protocol fi(t). 

In order to give a more definite interpretation of this 
quantity, let us consider, as in 1], the simple case of 
a particle bound to the origin by a spring of Hooke's 
constant k and to which a time-dependent force f(t) is 
applied. We then take / as a generalized coordinate, and 
write H(x, /) as 



H ( x > /) = \ kx2 - f x - 



Then 



and 



dH 
W 

dW = 



-x, 



-xdf. 



(14) 



(15) 



(16) 



On the other hand, the mechanical work that the external 
force exerts on the system is given by 



dW, 



mcch 



= fdx. 



Note that 



dW -dW mech = -d(fx). 



(17) 



(18) 



We see that the work appearing in Jarzynski's equality 
is not equal to the mechanical work exerted on the sys- 
tem. The difference between the two works is equal to 
the variation of the potential energy associated with the 
force. It is therefore no surprise that, if one evaluates 
the integral of the mechanical rather than the thermody- 
namical work, one obtains an incorrect estimate of the 
free-energy difference. 

We can now straightforwardly check the validity of the 
relation between the reversible work and the free energy, 
provided the full Hamiltonian is taken into account. We 
have, of course, 



dW mv = -df (x) 



(19) 



where (x)f is the equilibrium position of the particle at 
force /: 



(x) 



_!_ 

f~ k ' 



(20) 



Thus 



AG = -Jdf ( X } f = -Q = ~k{x) 2 f . (21) 

This quantity is equal to the one obtained from the par- 
tition function, as can be trivially checked. It is negative, 
because the positive variation of the spring potential en- 
ergy is more than compensated by the drop in the po- 
tential energy of the applied force. One can imagine, 
e.g., that the guide to which is constrained the particle 
is made rotated by an angle 8 around the origin in a ver- 
tical plane. Then / = mg sin 9, and, at equilibrium, the 
height of the particle has dropped to z — — (x)jSm0, 
with a corresponding gravitational potential energy drop 
equal to mgz. 

Of course nothing prevents us to take explicitly into ac- 
count this contribution, in order to evaluate the change of 
E — TS, where E is the internal energy. It is sufficient to 
subtract the potential of the applied force. This is indeed 
what is routinely made when estimating the free energy 
landscapes from manipulation experiments, following the 
suggestions of Hummer and Szabo ij] and several others. 

Let us now consider the second example discussed by 
Vilar and Rubf [lj , namely a harmonically bound particle 
subject to a harmonic time-dependent force: 



H(x,X) 



1 



~-kx 2 



1 



K{x - X) 



(22) 
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By expanding this expression, we obtain 

H(x,X) = ^{k + K)x 2 -KXx + \KX 2 . 
Then 



Of course 



dW = -dXK(x-X) 



(23) 



(24) 



In our case, (Ax 2 ) Q = (Ax 2 ) x . Thus 



(U)x ( U )o = 



l -K{{x) x -Xf 



Ik 2 



K 



K 



2 K 



(x) 



X 



(x) 



X 



(30) 



-dXK{(x) 



x 



X) 



(25) 



where {x) x is the average displacement of the particle 
when the center of the harmonic trap is placed at X: 



(x) 



K 



x 



K 



-.X. 



(26) 



Thus 

AG = 



dW rev = 



x 



dX' K 1 - 



K 



1 kK 



2k + K 



-X 2 



1 k(k + K) 

2 K 



k + K 

2 



X' 



x 



(27) 



This result is different from that reported by Vilar and 
Rubi, since it takes into account also the change in the 
potential energy U of the applied force. 

Let us evaluate the change of (U). We have 



(U)o = ^<* 2 ) = ^<A* 2 ) : 



(U) x = -KUx-X 



Ik 

2 



((x) x -X) 2 + (Ax 2 ) 



x 



(28) 



(29) 



We obtain therefore the difference AF of the internal 
free energy F = E - TS = ((H) - U) - TS: 



AF = AG — A (U) = ifc (x) 2 x . (31) 



We see that, with the correct application of the thermo- 
dynamic relations, it is possible to reconstruct the inter- 
nal (free) energy landscape. 

LP is grateful to Stan Leibler for pointing out to him 
ref. 0. 
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